Decay of Superfiow Confined in Thin Torus: 
A Realization of Tunneling Quantum Fields 



PL, 
On 



(N 
> 
O 



00 
OS 

C3 



c 

o 
o 



13 



M. Nishida and S. Kurihara 

Department of Physics, Waseda University, 3-Okubo, Shinjuku-ku, Tokyo 169-8555, Japan 

(February 1, 2008) 

The quantum nucleation of phase slips in neutral superfluids confined in a thin torus is investigated 
by means of the collective coordinate method. We have devised, with numerical justification, a 
certain collective coordinate to describe the quantum nucleation process of a phase slip. Considering 
the quantum fluctuation around the local minimum of the action, we calculate the effective mass 
of the phase slip. Due to the coherence of the condensate throughout the torus, the effective mass 
is proportional to the circumference L of the torus, and the decay rate has a strong exponential 
//-dependence. 



PACS numbers: 67.40.Hf, 64.60. Qb, 67.40.Fd 

The circulation of superfiow in a macroscopic torus 
container is quantized. However, it has been pointed out 
that the quantized circulation can decay via the thermal 
nucleation of phase slips in the case of very narrow torus, 
or extremely thin superconducting wires jjj-Q]. If the 
temperature is low enough, it is probable that the nucle- 
ation due to the macroscopic quantum tunneling become 
dominant ^|J^]. There have been several experiments of 
superconducting wires which, in fact, indicate the exis- 
tence of a crossover between the thermal activation be- 
havior in the high temperature regime and the quantum 
tunneling behavior in the low temperature regime . We 
expect that similar experiments on the neutral Bose gas 
systems will become available, in the near future, with 
the trapped atomic gases [||, or with the liquid 4 He in 
meso-porous materials ||. The purpose of this Letter is 
to estimate the probability of the quantum nucleation of 
phase slips in dilute Bose gas systems. 

At sufficiently low temperatures, a condensate in a di- 
lute Bose gas is described by a macroscopic wave func- 
tion tp(r). We consider a very thin torus container. The 
circumference L is much larger than the healing length 
£ = ; = where m is the mass of an individual atom, 
po is the density of the condensate (without flow), and 

_ inn a j g ^ interatomic interaction, with a being 
the s-wave scattering length. The cross-sectional area a 
is much smaller than £ 2 . We assume that the spatial de- 
pendence of ip in the cross section can be neglected by 
using the effective number density. We will discuss this 
point later. Taking x axis along the circumference of the 
container, the action S[ip] and the free energy F[ip] have 
the Gross-Pitaevskii (GP) form JHJ: 
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If we assume that the effective inertial mass of the 
phase slip is not greatly dependent on the tunneling path, 
the most probable tunneling path will be through the 
valley, passing through two local minima and one saddle 
point in the functional space of ip. The wave function ip 
is determined so as to minimize the free energy functional 
at each stage of the tunneling process. We calculated the 
minimum value of the GP free energy numerically with 
the constraint condition that the average velocity of the 
superfiow had each constant value. We adopt the units 
in which the unit of length is the healing length £, the 
unit of velocity is the sound velocity c = y/^-, and the 
unit of number density is the density of the condensate 
Po- A dimcnsionless parameter, a = Po<t£,, determines 
the importance of quantum fluctuations |5| . We assume 
a ^> 1, and use WKB approximation. Considering a <C 

£ 2 , a -C po£ 3 oc Pg 1 ^ 2 , so a is large when the density is 
low. We include the following term into the free energy 
Eq. (|): 
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where p(x) and 9(x) are the density and the phase of the 
condensate respectively, v s is the constraint value of the 
superfiow velocity, and A is supposed to take sufficiently 
large value which is determined by the condition that the 
difference between the constraint value v s and the aver- 
age velocity of the obtained solution be smaller than, say, 
1.0 x l(T 6 c. We had L = 20 £ in a mesh of 100 points, 
and used the simulated annealing method pl| . Because 
there are a lot of local minima which correspond to mul- 
tiple phase slip solutions, it is difficult to obtain true 
minimum values with other methods. In order to avoid 
being trapped in these local minima, we introduce ran- 
dom noises which correspond to thermal fluctuations at 
given temperature. The calculation proceeds gradually 
from high temperature region to low temperature region, 
and continues until equilibrium values are reached. In 



our calculation, the final temperature is 1.5 x 10 



-10 



where p ~ pog is the chemical potential. 



Po9, 



and the expected error of the free energy is O (10 4 pog) ■ 
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Extrema of the GP free energy functional are solutions 
to the time-independent GP equation: 
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There are two kinds of solutions to this equation || . The 
first kind of solution corresponds to a free-energy local 



minimum: 
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where k is the phase gradient or, in other words, the 
superflow velocity. The second kind corresponds to a 
free-energy saddle point: 
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The phase gradient of if) B approaches k asymptotically as 
we go away from the phase slip center {x = 0). Here- 
after, we call k "asymptotic velocity" . If we can find the 
function which interpolates these two kinds of solutions 
smoothly, and fits the numerical results, we may use it as 
the wave function which describes the condensate during 
the whole nucleation process. Considering the periodic 
boundary condition, we have found the following inter- 
polation function: 



x < sin - tanh 
I 2 



— ain-x +icos- . (7) 



We must determine the parameter q so that ipi satisfies 
the periodic boundary condition. If qL 3> 1, we can take 
the form, 



q = kL — 27rn, 



(8) 



where n is the winding number of the phase around the 
circumference of the container (n = 0, ±1, ±2, . . . ). The 
parameter q expresses the phase difference between two 
sides of the phase slip center. Hereafter, we call it "phase 
jump" . We can adopt either asymptotic velocity or phase 
jump as an appropriate collective coordinate to describe 
the nucleation process. Substituting Eq. ^ into Eq. (||), 
the potential energy for the asymptotic velocity k is ob- 
tained as follows: 
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Using fitting parameter, we fit the amplitude and 

phase of the wave function Eq. M) with those which are 



given by the numerical minimization. The phase jump 
changes about L(= 20)% when k changes 1% (see Eq. 
(||)), so the shape of the phase is very sensitive to the 
value of k. Fig. |l| shows that Eq. (||) is in good agreement 
with the numerical results. 
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FIG. 1. Potential energy for the collective coordinate k. 
The points are the minimum values of the GP free energy 
which were calculated numerically. The solid line denotes the 
value of the potential energy Eq. (^) which is calculated using 
the interpolation wave function Eq. ([?]). 



We then consider the effective inertial mass for the mo- 
tion of the collective coordinate. Because of the trans- 
lational invariance, the total momentum must be con- 
served between the initial state and the final state. If 
there were no contributions from the noncondensed com- 
ponent, the condensate would not lose the momentum 
and stay at the local minimum. However, quantum fluc- 
tuations of tp exist around the distorted condensate wave 
function ipi, and these fluctuations will play the role of 
a noncondensed component. Therefore, the condensate 
can oscillate around the local minimum and tunnel to an- 
other local minimum by exchanging the momentum with 
quantum fluctuations, even at zero temperature. The fre- 
quency of the zero-point oscillation of the collective co- 
ordinates around the local minimum will have the same 
value of the collective mode whose wave length is L. Be- 
cause we have already determined the potential energy, 
Eq. (^), we can estimate the effective inertial mass from 
the frequency of the collective mode and the curvature of 
the potential energy for n or q at the local minimum. In 
this procedure, we assume that the inertial mass is not 
greatly dependent on the collective coordinate, and ap- 
proximate it to the value at local minimum. The coupling 
with the collective modes also suggests that there are the 
effects of dissipation, which will suppress the tunneling 
probability Jl2] |. However, we will not take into account 
these effects, but concentrate our interests on the bounce 
motion of n or q. In order to clarify these two points, we 
would have to consider quantum fluctuations explicitly 
and integrate them out, thus requiring further study. 

In order to estimate the frequency of the collective 
mode, we consider the fluctuation around the local min- 
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imum solution, Eq. (||). We substitute ip = ip m + (u + 
iv)e 1KX into the action, Eq. (|j), and retain the terms up 
to the second order of u, v, then we have 

S[u,v]=Ha Jdt -V(k) + J dxi.2iv {^sj 1 - -7Q 



(u v)G~ 1 



(10) 



where G is the Green's function for u, v, whose Fourier 
component is given by 
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G has the poles which correspond to the Bogoliubov 
modes in the flowing medium. We can derive these poles 

by Galilei transformation from the co-flowing frame to 

2 

the lab frame, noting that the density is now 1 1-. 

From these poles, we estimate the effective angular fre- 
quency of the oscillation of k by the following relation: 



eS 



■ UJ 



(12) 



rate F is given by T — Ae~ SB / h , where Sb is the bounce 
action, which is the Euclidean action for the bounce mo- 
tion of k [Ol : 



ha / dr 



M K (n) . 



k 2 + V(k) - V( Kn ) 
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Because k contains either the winding number n or the 
system length L, it is more convenient to use the phase 
jump q for the collective coordinate when we consider 
the dependence of the tunneling rate T on the superflow 
velocity n n or the system length L. Rewriting the bounce 
action Sb by using the phase jump q = (k — K n -i)L, we 
obtain 



2ha 



dq^2M q (n)V n (q), (17) 



where qo is the classical turning point of — V n (q), 
M K (n) L 



M q {n) 



L 2 



(L » 1), (18) 



V n {q)=V(q/L + K n - 1 )-V(K n ) 
,2 
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Around the nth local minimum, whose velocity is K r , 
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j-, the effective angular frequency becomes 



1 - -< 



L 
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When q is so small that we cannot assume qL ^> 1 , Eq. 
@ and Eq. (^) become inaccurate and we must deter- 
mine the relationship between q and k with direct consid- 
eration of the periodic boundary condition. After doing 
so, the potential curvature at the nth local minimum be- 
comes 



V"{K n ) = L [ 1 - -K 2 n ) , 



(14) 



where the prime denotes differentiation with respect to 



From Eq. (13) and Eq. (|14|), the effective inertial mass 
for k is given by 



M K {n) = — - 



1 - ^K 2 
1 2 K n 



L 3 



(2ny 



(L » 1). 
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We now consider the tunneling problem from nth local 
minimum to the (n— l)th local minimum. The tunneling 



We thus find that the effective inertial mass for the nu- 
cleation of a phase slip is proportional to the system 
length L. It basically comes from the fact that the po- 
tential curvature is proportional to the system length, 
(see Eq. (|l4|)). When a phase slip nucleates, the whole 
wave function of the condensate changes its configuration 
due to its overall coherence, so the free energy change be- 
comes larger as the system circumference becomes longer. 
Therefore the coherence of the condensate wave function 
is the origin of the L-dependence of the inertial mass. 

We can calculate the bounce action analytically in the 
zero-current limit (k„ — > 0): 



_ 2V2haVL f 2 * , .2 
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Fig. U shows the bounce actions normalized by So 
in the five different cases with the circumferences of 
50 — 10000 £. In all cases, the ratio Sb/So goes to zero 
at the common velocity v c ~ y/2/3c. This velocity cor- 
responds to the critical velocity of the superflow. If the 
asymptotic velocity k goes beyond v c , the saddle point 
solution to the time- independent GP equation, Eq. (|^), 
will vanish. No stable local minima exist for velocity 
n n > v c . The value of v c agrees with the one which is de- 
rived by Landau's criterion. As K n approaches the value 
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FIG. 2. Bounce actions normalized by So, which are calcu- 
lated for the five different samples with the circumferences of 
50 - 10000 £. 



of v c , the bounce action becomes extremely small. In the 
case of L = 100£,n = 12, the normalized bounce action 
Sb/S becomes 0.012. 

We see from Fig. || that there is no sizable L- 
dependence for L>1 under consideration. If L is larger 
than 100 £, all the data seem to be on the same line. 
Therefore we conclude that the L-dependence mainly 
comes from the effective inertial mass, and the bounce 
action is proportional to L 2 . This situation is completely 
different from that in the thermal nucleation processes. 
In the thermal nucleation processes, the decay rate is es- 
timated by the Arrhenius formula T oc e - SF / k ^ T ^ which 
is determined only by the barrier height SF, so there is 
no dependence of the decay rate on the system length 
In the quantum nucleation, the coherence of the 
condensate wave function becomes important, and the 
decay rate has the strong L-dependence: 



T oc exp{-1.82a\/Z/(K„)}, 



(21) 



where f(K n ) is the function which depends only on the 
superflow velocity, and whose shape is shown in Fig. ^. 

We now discuss the spatial dependence of ip in the 
cross section. If we consider a complete set of radial 
eigenfunctions, the energy scale of the ^th eigenstate will 

fc2 2 *-2 

become . Because it is much larger than u ~ in 

ma ° ' rat/ 1 

our model when v ^ 0, the v = channel will become 
dominant [Q. Besides, the energy scale of the potential 
barrier is a/i ~ , which corresponds to the energy 

2 

of 7^3- particles excited into the vth channel. The ratio 



to the total number is 



([HP 



< 1 when v ^ 0. These 



estimations seem to indicate that we can safely neglect 
the spatial dependence in the cross section. However, 
there are other effects to be considered: For example, 
the nonlinearity of GP equation. It will cause couplings 
between the current direction and the radial direction. 
If we take this into account, we may find a more appro- 
priate tunneling path. Moreover, when the condensate 
is squeezed in a region whose size is smaller than the 



healing length, the phase of the condensate will fluctuate 
intensively. Because phase fluctuations (((56*) 2 ) will ef- 
fectively suppress the amplitude of the wave function by 
the factor 1 - \ {(S9) 2 ) +O ((S9) 4 ). This means that the 
effective number density of the condensate, and thus the 
parameter a, which gives the scale of bounce action, will 
become smaller than the results of the present calcula- 
tion. Therefore, it is probable that our result of tunneling 
rate is an underestimate. The analysis of these effects is 
left for the future work. 

In conclusion, we have studied the quantum nucleation 
of phase slips in torus shaped superfluids by means of 
the collective coordinate method. We have found a wave 
function which interpolates the local minimum and sad- 
dle point solutions to the time-independent GP equation, 
and this fits our numerical results. We have determined 
the appropriate collective coordinate which describes the 
whole nucleation process, and obtained the analytic form 
of the potential barrier. We have shown that the effective 
inertial mass is proportional to the system length L due 
to the coherence of the condensate. Because the length 
dependence of the bounce action mainly comes from the 
mass, we have concluded that the bounce action is pro- 
portional tO Zj2. 
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